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We study the behaviour of the inverse participation ratio and the localization transition in in- 
finitely large random matrices through the cavity method. Results are shown for two ensembles of 
random matrices: Laplacian matrices on sparse random graphs and fully-connected Levy matrices. 
We derive a critical line separating localized from extended states in the case of Levy matrices. 
Comparison between theoretical results and diagonalization of finite random matrices is shown. 

PACS numbers: 



I. INTRODUCTION 

Since the pioneering work of Wigner in nuclear physics 
[l| , random matrices have been intensively studied due to 
their wide range of applications in several fields of physics 
and other disciplines. Some examples include quantum 
chaos [31, localization in electronic systems 3, diffusion 
in random graphs [J, la] , finance [g, Q and complex net- 
works Q. 

In random matrix theory, one is interested in physical 
quantities that can be computed from the eigensolutions 
of a sample drawn from an ensemble of TV x TV random 
matrices. One of the main quantities of interest is the 
density of states (DOS). In the case of the Gaussian or- 
tliogonal ensemble the DOS obeys the Wigner semicircle 
law in the limit iV — > oo [9| . There are several ensembles 
in which the average DOS differs from the semicircle law 
in a nontrivial way. The examples include the ensemble 
of sparse matrices ll(H 15 | . Laplacian matrices [3, Ila-[l3 
and Levy matrices pOl - [23 | . Finite size effects do not play 
a significant role and the DOS converges to its large N 
limit relatively fast. 

Another important quantity is the inverse participa- 
tion ratio (IPR) since it provides valuable information 
about the nature of the cigcnstatcs. The IPR allows one 
to quantify the number of nonzero components in a cer- 
tain eigenvector in the limit N ^ oo. It is a suitable 
parameter to describe quantitatively a delocalization- 
localization transition, since it distinguishes between 
eigenstates that have a finite number of nonzero com- 
ponents (localized states) and eigenvectors that have 
an extensive number of nonzero components (extended 
states). The critical eigenvalue separating localized from 
extended states is called the mobility edge or the localiza- 
tion threshold. The mobility edge determines the Ander- 
son transition in electronic systems [3| and the emergence 
of particle traps in diffusion models on random lattices 
[J, Q . An equation that determines the mobility edge in 
the ensemble of sparse random matrices was obtained by 
means of the supersymmetric method 23[. Within our 



knowledge there exists no full numerical solution for the 
mobility edge by calculating the IPR with the supersym- 
metric method. 

In contrast to the DOS, numerical diagonalization re- 
sults show a significant dependence of the average IPR 



upon TV [ij, [l^, [l3l ■ Moreover, since localized states are 
usually present in the tails of the spectrum of random 
matrices, one would have to diagonalize extremely large 
matrices to detect these states. In order to determine the 
localization threshold in the limit A^ — > oo, the imaginary 
part of the self-energy Q and the variance of the den- 
sity of states [21, [23] have been proposed as appropriate 
parameters. However, a more quantitative description of 
the localization transition would be obtained by calculat- 
ing the average IPR for infinitely large random matrices. 

Laplacian matrices on sparse random graphs |16| and 
fully-connected Levy matrices [20] are two examples of 
random matrices in which results for the average IPR 
as a function of the eigenvalue have been obtained only 
through diagonalization of finite matrices. Laplacian ran- 
dom matrices arise, e. g., in the study of diffusion on 
random graphs [j, [^, \1^ and in the instantaneous normal 
modes approach for liquid dynamics [17{ . Levy random 
matrices appear, e. g., in models of spin-glasses with 
dipolar RKKY interactions [2g, in the study of disor- 
dered electronic systems with interactions decaying as a 
power-law of the distance |2C| , in portfolio optimization 
[271 ] and in the study of correlations in data, for instance, 
coming from financial time scries [2a [23] ■ 

In this work we calculate the eigenvalue-dependent 
IPR from the Green function corresponding to the ran- 
dom matrix in the limit A^ — > 00. The cavity method 
provides a self-consistent equation for the Green function, 
which is easily solved through a population dynamics al- 
gorithm [sol- This approach is free of finite-size effects 
and it can be used to study the average IPR of differ- 
ent ensembles of infinitely large random matrices. We 
present results for the average IPR as a function of the 
eigenvalue and the presence of a localization transition 
for Laplacian matrices on a sparse random graph and for 
fully-connected Levy matrices. In the case of Levy ma- 
trices, we calculate a critical line that separates localized 
from extended eigenstates. The theoretical results for 
N ~^ 00 arc compared with numerical diagonalization of 
finite matrices and with results of previous works. 

In the next section we define the average quantities 
of interest and we discuss how they can be calculated 
in the limit A^ — )■ 00. The results for the average IPR 
and the localization transition in the case of Laplacian 
and Levy matrices are shown in section IIIII In section 



IIVI wc present our conclusions. An equation that relates 
the Green function and the IPR is derived in appendix 
lAl We include a detailed discussion of the cavity method 
for the ensemble of Levy matrices in appendix [BJ 



II. THE GENERAL SETTING 

In this section we show how the DOS and the IPR are 
written in terms of the diagonal elements of the Green 
function. The joint distribution of the real and imaginary 
parts of these diagonal elements is the central quantity 
of interest to determine the average DOS and IPR. We 
discuss how the cavity method can be employed in order 
to calculate this distribution in the large N limit. 



which is the limit iV — >■ oo of the average value of Y^ 
over all the states lying between A and X + dX when p{X) 
is finite. 

The quantity that allows one to determine the DOS 
and the IPR is the Green function associated to J. Its 
diagonal elements are defined as follows 
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where z is the complex variable z = X — ie. The DOS can 
be obtained in the limit e — > 0+ according to 
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A. Random matrix parameters 

We define an ensemble fl^ of symmetric N x N ran- 
dom matrices with real elements. Assuming that a given 
matrix J G ^n has a set of eigenvalues {A^}^i=i....^7v and 
normalized eigenvectors {| /i)}^i=i,...,Ar, the DOS lying 
between A and A + dA is given by 
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fl = l 



The IPR associated to a given eigenvector | fi) is de- 
fined as 
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where V'u ^ (* I p) is the component i of the eigenvector 
I p). The set of normalized vectors {| i)}i=i,...,Ar is the 
canonical site basis. The IPR allows one to distinguish 
between two extreme situations in the limit iV — >■ oo. In a 
delocalized or extended region of the spectrum, a number 
of sites of 0{N) contributes to a given eigenstate. Due 
to the normalization X]i=i(''/'u)^ = 1' the components 
of a given | y) satisfy -0^ = 0{\/^/N) and, as a conse- 
quence, we have limTv^oo Y^ = 0- In contrast, only a 
finite number of components {i/'u} is nonzero for a local- 
ized eigenstate. The components of a state /x localized 
on d sites satisfy ip^ ~ 0{l/^/d) and the IPR is given by 
Y^ ~ 0{l/d) in the limit A^ — > oo. This simple analysis 
shows that the IPR is a suitable parameter for a quanti- 
tative description of the transition between extended and 
localized eigenstates. 

Here we are interested in the average behaviour of Y^ 
over all the states in an infinitesimal region of the spec- 
trum. We define the eigenvalue dependent IPR 



The IPR is expressed by 



P{X) == lim lim -- 
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in the regions of the spectrum where there are no de- 
generate states. Equation © has been employed in the 
study of localization properties through the supersym- 
metric approach |23j . In appendix [^ we explain how to 
derive eq. ([6]). Finally, we rewrite eqs. ([5]) and ([6]) as 
follows 
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where we have introduced the average (/(w)) = 
/ da;WA,e(w)/(w) with respect to the joint distribution 
Wa,£(w) of the real and imaginary parts of G^{z) 
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We have assumed that the distribution of G^{z) is a 
self-averaging quantity in the limit N -^ oo, with {. . .) j 
denoting the ensemble average with respect to the distri- 
bution of J. This implies that p(A) and P{X) are self- 
averaging quantities in the limit N ^ oo. The integrals 
present in (/(w)) run over the entire real and imaginary 
axes of uj. Our goal consists in calculating the distribu- 
tion Wa,c('^), since this quantity allows us to compute 
the average DOS and IPR through eqs. ^ and ([5]). 



A cavity calculation of the distribution of Green 
functions 
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In this subsection we explain how the cavity method 
(3) can be used in order to calculate WA,e(w). As an exam- 
ple, the technical details involved in the calculations are 



shown in appendix [B] for the ease of fuUy-connected Levy 
matrices. 

The diagonal elements of the Green function can be 
written as a Gaussian integral over the real vectors x = 
{xi,...,xn), namely 
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Defining the normalized complex function 
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with 
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the elements G^j^{z) assume the form 



Gkki^) = i / dxxlVN,z{x) = i / dxkxlVN,z{xk), 

(13) 
where VN,zixk) follows from Vn,z{x) by integrating over 
all variables besides Xk- More generally, we write for a 
set A of indices 
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of the variables in the neighbourhood di factorizes on the 
cavity graph G'-'-' , i.e. 



N,z 



(Xd,) 



n 



N,z 



(^,) • 



(15) 



The cavity graph G'*-' is the subgraph of G where node 
i and all its connections have been removed. The 
marginals {'P^\,{xi)} are defined on the cavity graph. 
Since the local marginals {'Pn,z{xj)} on the real graph 
can be written in terms of the local marginals on tj('), 
one can solve first the problem on the cavity graph as 
a function of {r^ zi^j)} ^^d then reconstruct the local 
marginals {Vn,z{xj)}. 

In the case of disordered spin systems defined on 
finitely connected graphs, the cavity method is also 
known as the Bethe-Peierls iterative method [301 ■ ^^ Ihi^ 
case, we expect that the factorization assumption p5)) 
holds outside the spin-glass phase since the graph looks 
locally like a tree. When the graph is a tree the condi- 
tion P^ holds for iV — ^ oo. In the case of disordered 
spin systems defined on fully-connected graphs, the van- 
ishing of the connected correlation functions for N -^ oo 
ensures that the assumption (fT5|) holds outside the spin- 
glass phase [3l| . 

By employing the cavity method and following analo- 
gous calculations as done in [l3|, we show in appendix [B] 
that the marginal VN.z{xk) at site k is given by 
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Equation (fT3|) shows that the Green functions {G^(z)} 
can be computed for a single instance of J once we know 
how to calculate the local marginals {'PN,z{xk)}- 

Disordered spin models on a random graph are defined 
through a probability measure on this graph. The cavity 
method provides an efficient way to compute the local 
marginals of disordered spin systems defined on fully- 
connected [3 if and finitely connected random graphs [30[ . 
The application of the cavity method to the study of ran- 
dom matrices relies on the analogy between the random 
matrix problem and disordered spin systems defined on 
graphs [iJl ■ Using this similarity one can apply the cav- 
ity method to calculate the marginals {VN,z{xk)}- 

One can associate a random graph to a random matrix 
as follows: the graph contains N nodes. A certain vari- 
able Xi {i = 1,. . . ,N) is associated to the corresponding 
node i of the graph. When J^ = there is no edge be- 
tween nodes i and j, while this pair of nodes is connected 
when Jij 7^ 0. The interaction strength between nodes i 
and j is given by the value of Jij, such that J specifics ^j^g form 
the topology of the random graph and the interaction 
strengths between the nodes. 

For a given graph instance G, di is the the set of nodes 
connected to a certain node i. The cavity method is 
based on the assumption that the distribution 'Pn,z{xq.) 



The diagonal elements of the Green function are deter- 
mined from the fixed-point solution of the equations 
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for i = 1, . . . ,N and for all k £ di, where di is the set 
of all indices j in a given row i such that Jy ^ 0. The 
symbol di \ k denotes the set di without index k. The 

quantities {G^^'^ (z)} are the diagonal elements of the 
Green function of the matrix following from the original 
matrix through removal of row k and column k. The 
specificity of the random matrices under study can man- 
ifest itself only in the number of indices present in di and 
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in the form of the function hfJJij, G^.- 

For fully-connected Levy matrices and Laplacian ma- 
trices considered in this work, h^- (Jij , G,,, (z)) 
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Levy matrices: 



hij[Jij,G^j (2)) — JijGjj (2) , 



N,iz) 



(19) 



Laplacian matrices: 



,Af, 



T r' 



N.ii 



(^)) 



JijGjj (z) 



1 + J^■,G 



JV,W/ 



J, 



tj ■ 



(20) 



Equations (HZl) and dUl) with hfj{Jij,G]f'{z)) given by 
P^ have been obtained previously in the study of fully- 
connected Levy matrices [20[ and sparse random matrices 

A self-consistent equation for WA,e(w) is obtained by 
substituting eq. (|18p in eq. © and performing the av- 
erage over the ensemble of random matrices. We have 
solved numerically this self-consistent equation through 
a population dynamics algorithm |30[, which consists in 
parametrizing the distribution Wx^ei'^) by a large popu- 
lation of stochastic variables representing instances of w. 
At each iteration step, one of these variables is chosen at 
random and updated according to its probability distri- 
bution, until a stationary form for Wx^e{uj) is reached. A 
detailed discussion of the population dynamics method 
in the context of random matrices and the corresponding 
algorithm are presented in 
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According to eqs. ([7]) and (|H]), one has to obtain nu- 
merical results for the distribution W\,e{i^) in the limit 
e — >■ 0. One has to calculate WA,e(w) for very small but 
finite values of e, since Dirac delta peaks might arise, 
for instance, in the spectrum of sparse random matrices. 
In this way the Dirac delta peaks are approximated by 
Lorentzian functions with a finite width e [13l [ij] . In the 
next section, we specify the ensembles of random matri- 
ces and the corresponding distribution W\^t{uj) for each 
case. 



III. RESULTS 

In this section we show the results for two different 
ensembles of symmetric random matrices: Laplacian ma- 
trices on sparse random graphs and fully-connected Levy 
matrices. In both cases we focus on the behaviour of the 
average IPR and the presence of a localization transition. 



with Cii = for Vi. In the limit N -^ 00, J has a sparse 
structure and the number of nonzero elements per row 
exhibits a Poissonian distribution with average c. The 
nonzero elements {Kij} are drawn according to the dis- 
tribution pK{Kij). We consider here two different cases: 
(i) the elements {Kij} assume a fixed value Kij ~ — 1/c, 
such that pK (Kij) = S{Kij + l/c); (ii) the elements {Kij} 
arc drawn from a Gaussian distribution with zero mean 
and variance 1/c. 

Inserting eq. (jlSp in eq. Q and performing the ensem- 
ble average, one can derive the following self-consistent 
equation for W\,e(w) 
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The updating function H{uj,K) is given by 



H{lu,K) = 
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1 + Kuj 
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(23) 



(24) 



The population dynamics algorithm can be used to solve 
eq. ([23]) numerically |30l |. 
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A. Laplacian matrices 

The elements of the Laplacian matrix J^ on a random 
graph can be defined according to jl4 



N 



^ij — ^ij^ij ^ij / J ^ik-^ik 



(21) 



fc=i 



in which c and K are symmetric matrices. We consider 
here only the case in which the elements of the connec- 
tivity matrix c are i.i.d.r.v drawn from the distribution 



Pc{Ctj) 



l-^)'5c.,o + ^^c.„i, (22) 



FIG. 1: Population dynamics results for the average DOS 
(solid lines) and IPR (dotted lines) for Laplacian matrices 
with c — 20 and elements {Ki} drawn from a Gaussian dis- 
tribution. These results were obtained with e = 0.001 and a 
population of 5 x 10^ samples for the distribution W\^t{i-^)- 
Numerical diagonalization results for the DOS (*) and the 
IPR (o) obtained with an ensemble of 500 matrices of di- 
mension A'^ = 3000 are shown. Error bars for the IPR are 
indicated. The insets show population dynamics results for 
the DOS (solid lines) and the IPR (dotted lines) as a function 
of A in the left and right tails of the spectrum of Laplacian 
matrices with c = 20, e = 0.001 and fixed elements {Ki}. 

In fig. [T]we illustrate the results for p(A) and P{X) ob- 
tained from the population dynamics algorithm and from 



diagonalization of finite matrices for c = 20. The spec- 
tra of random matrices with a sparse structure contain 
delta peaks located at the eigenvalues of isolated finite- 
size clusters for any value of c [ll|, ll3, 13 • Eq. ^ gives 
an approximation for the IPR in these regions of the spec- 
trum due to the presence of degenerate states. In order 
to minimize the effect of these singular contributions, we 
have chosen a large value of c. The main graph of fig. [1] 
illustrates p{X) and P(A) for the Gaussian distribution of 
{Ki} only when A > 0, since the spectrum is symmetric 
around zero. The insets show the behaviour of p(A) and 
P(A) in the tails of the spectrum for the case of fixed 
Ki = — 1/c. In both cases, P{X) is vanishingly small in 
the central part of the spectrum, corresponding to a re- 
gion of extended eigenstates. The eigenvectors undergo 
a localization transition in the tails of the spectrum, as 
shown by the increase of P(A). It is difficult to determine 
the IPR through numerical diagonalization in the tails of 
the spectrum since one has to diagonalize extremely large 
matrices. 

By means of numerical diagonalization and a siiigle 
defect approximation (SDA) [i37llali the authors of [ia| 
have studied the eigenstates corresponding to the regular 
peaks that appear for large and small eigenvalues in the 
spectrum of Laplacian matrices with fixed values of {Ki}. 
They have found that these states are localized on a finite 
number of sites that have a small or large connectivity 
in comparison to the mean c. Fig. [T] complements these 
results by showing that P(A) -^ 1 for large values of 
A, which means that eigenstates corresponding to large 
eigenvalues are localized on a single site. The presence 
of the peaks is reduced when one introduces Gaussian 
disorder in the elements {Ki}. 

The agreement between diagonalization and theoreti- 
cal results for p(A) in fig. [T]is very good. We have found 
that p(A) depends weakly on N or e in the case of numeri- 
cal diagonalization or population dynamics, respectively. 
In the case of f (A), both results exhibit a very good 
agreement in the central region of the spectrum and in 
parts of the tails, for the particular values of e and TV 
chosen. The results show a discrepancy in the far regions 
of the tails where the eigenvalues are rare, as shown by 
the higher fluctuations in the numerical diagonalization, 
illustrated by the error bars. However, the average IPR 
depends on the values of N and e. 

Figure [5] illustrates the behaviour of P(A) as a func- 
tion of N and e for the Laplacian matrix with Gaussian 
elements {Ki} and c = 20. The results show that, for 
A = 2.90 and A = 3.35, the average IPR goes to zero 
when e — >■ 0. Accordingly, the diagonalization results 
for P(X) exhibit a similar qualitative behaviour as N in- 
creases. For A = 6.70, the average IPR has a finite value 
for e — )- 0, since this is the region of localized eigenstates. 
These results indicate that the localization transition pre- 
sented in fig. [1] becomes sharper for N ^ oo. By em- 
ploying numerical diagonalization methods, the authors 
of [l6| have found the value Ac — 1.67 for the mobility 
edge on the right tail in the case of Laplacian matrices 
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FIG. 2: Average IPR for Laplacian matrices as a function of 
e (population dynamics, top graph) and A'^ (numerical diago- 
nalization, bottom graph) for c = 20 and different numerical 
values of A, which are shown explicitly on the figure. The 
elements {Ki} are drawn from the Gaussian distribution de- 
scribed in the text. The diagonalization results were obtained 
considering an ensemble of 500 matrices for each A'^, with error 
bars included in this case. The population dynamics results 
were obtained with a population of 5 x 10^ samples. 



with c = 20 and fixed elements {Ki}. We have calculated 
approximately the localization threshold Ac on the right 
tail by using the e independence of P{X) as a criterion to 
determine the localized region. For c = 20 we have found 
the values Ac — 1.95 and Ac — 5.65 for Laplacian matrices 
with fixed and Gaussian elements {Ki}, respectively. 



B. Levy matrices 

The fully-connected Levy matrix is a symmetric matrix 
in which Ja ~ for Vi. The nondiagonal elements are 
i.i.d.r.v. drawn from the Levy distribution Pa{J), defined 
through the characteristic function La{q) 
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The characteristic function is of the form 
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The distribution Pa{J) is fully determined by the param- 
eter a S (0, 2]. For a < 2, a characterizes the power-law 



decay of Pq. ( J). We consider only Levy distributions cen- 
tered around zero. The scaling with N in eq. (PS|) ensures 
that the spectrum converges to a stable form in the limit 
TV — >■ oo |20l- The distribution of Green functions does 
not depend on the skewness parameter [20|, ^^ ■ 

For a ~ 2 wc recover the Gaussian orthogonal ensem- 
ble since Pa (J) is a Gaussian distribution with zero mean 
and variance 1/A^. For a < 2, the asymptotic behaviour 
of Pa{J) for I J| — )■ oo can be derived from the explicit 
form of La (q) '■ 
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The integrals for the second and higher moments of the 
distribution diverge for a < 2 due to the power-law decay 
iUustratcd by eq. ([ST]) . 

Due to the power-law tails of the Levy distribution, 
each row of J contains an infinite number of elements 
of order 0{N~^^°') and a finite number of elements of 
order 0(1). For small values of a, it has been argued 
that fully-connected Levy matrices can be seen as sparse 
random matrices [2Cl |21| . The spectrum of Levy matri- 
ces has been calculated with the cavity method in pre- 
vious works |20l . |22| and an equation for the distribution 
of Green functions has been determined using the gen- 
eralized central limit theorem [32] • We have followed a 
different approach to calculate the distribution of Green 
functions, in which the underlying sparse character of 
Levy matrices becomes transparent. Besides that, the 
resulting self-consistent equation can be solved through 
a population dynamics algorithm, which is a practical 
advantage in comparison with previous works where one 
has to deal with a complicated system of integral equa- 
tions [22|. 

In order to take the ensemble average and the limit 
A^ — >■ (X) of the distribution of Green functions, we intro- 
duce a cutoff 7 that makes an explicit distinction between 
strong matrix elements Jij > 7 and weak matrix elements 
Jij < 7 (see appen dix iBt. This trick has been introduced 
in spin systems in |33l . |3J] The backbone of strong ma- 
trix elements can be treated as a sparse random matrix, 
leading to the following self-consistent equation 



H{uj, J) = J'^uj 



PJniJ) 



and 



2|J|° + 1 





2Ca 



|J|>7 
|J|<7 



a7" ' 
27'-"Ca 



(30) 

(31) 

(32) 
(33) 



The distribution of Green functions follows from 
M^A,c(w) = lim^^o W^A,<r,7(w). 

The quantity c-y is the average number of strong matrix 
elements per row and p,j^-y{J) denotes their distribution. 
The contribution of the infinite number of weak matrix 
elements is taken into account through the law of large 
numbers, leading to a term proportional to their vari- 
ance ai. Eq. (|29p shows that the distribution W\_^_^{uj) 
contains a part coming from a sparse random matrix of 
strong matrix elements and an average contribution due 
to the weak matrix elements. 

We have solved eq. (PH]) through a population dynam- 
ics algorithm. The idea is to obtain results for small 
values of the cutoff 7. In fig. [3] we show results for the 
DOS of Levy matrices obtained from the numerical di- 
agonalization of finite matrices and from the numerical 
solution of eq. ((29)) . The DOS of Levy matrices is sym- 
metric around zero. For both values of a the agreement 
between diagonalization and population dynamics results 
is excellent. 
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FIG. 3: Comparison between numerical diagonalization (full 
lines) and population dynamics results (different types of sym- 
bols) for the average DOS of Levy matrices with a = 0.75 
(*) and a — 1.25 (o). The diagonalization results were ob- 
tained considering an ensemble of 1000 matrices of dimension 
N — 1500. The population dynamics results were obtained 



considering e ; 
pies. 



0.001, 7 = 0.01 and a population of 10 sam- 
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FIG. 4: Population dynamics results for the average IPR of 
Levy matrices as a function of 7 for a = 0.5, A = 1 and a 
population of 5 x 10^ samples. Three different values of e are 
shown: e = 0.01 (o), e = 0.001 (x) and e = 0.0001 (*). The 
average IPR depends on e for 7 — >■ 0. 
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FIG. 5: Population dynamics results for the average IPR of 
Levy matrices as a function of 7 for a = 0.5, A = 5 and a 
population of 5 x 10® samples. Three different values of e are 
shown: e = 0.01 (o), e = 0.001 (x) and e = 0.0001 (*). The 
average IPR is independent of e for 7 — >■ 0. 



For a — >• 2, we obtain c^ -> and cr^ — > 1, and the 
distribution Wx^e,-y{'^) reduces to the simple form 



WxA^)^s[^- 



1 



(34) 



When inserted in the definition of {to) , eq. p4p gives rise 
to a quadratic equation for (w) vifhose solution leads to 
the Wigner semicircle law by means of eq. d?]) . 

In figs, m and [5] we illustrate the population dynamics 
results for the behaviour of ^(7) as a function of 7 for 
A = 1 and A = 5, respectively. In both figures we consider 
a — 0.5 and three different values of e. Figure |4] shows 
that for 7 — > the average IPR decreases for decreasing 
values of e, which is an indication that A = 1 corresponds 
to a region with delocalizcd states. As 7 decreases in fig. 
[5l -P(7) converges to a finite value that does not depend 
on e, which corresponds to a region of localized states. 

We have used the e independence of -P(A) in the local- 
ized region as a criterion to calculate approximately the 
localization threshold in the limit N — )■ 00. In fig. IHlwe 
present results for the critical line separating localized 
from extended states in the (a. A) plane for 7 = 0.008. 
For small values of a, c-y is small and the tails of the dis- 
tribution pj^j{J) are very long. In this case the sparse 
matrix character of Levy matrices is highlighted and the 
region of localized eigenstates is larger. For increasing 
values of a, the parameter c^ increases and the tails of 
the distribution pj.-yiJ) are less long. The distinction 
between strong and weak matrix elements becomes less 
important and the fully-connected character of Levy ma- 
trices is highlighted, leading to a larger region of extended 
eigenstates. The results for the localization transition 
when a is large are very noisy due to the larger values of 
Xr involved in the calculations. This makes the results 



for a > 1.3 very inaccurate. Besides that, the population 
dynamics algorithm becomes slower for increasing values 
of Cy. We have obtained numerically that the average 
number of strong matrix elements c-y reaches its maxi- 
mum value at a ~ 2 when 7 — >■ 0. This indicates that 
the region of extended states is the largest possible for 
a~2. 

For a = 0.5 the population dynamics results show that 
P(A) -^ 1/2 as A ^ 00. This result agrees with the dis- 
cussion presented in [20|. According to this work, due 
to the strong fluctuations of the Levy matrix elements, 
the eigenstates corresponding to large eigenvalues are lo- 
calized on pairs of very strongly interacting sites, which 
leads to an IPR equal to 1/2. Based mostly on numerical 
diagonalization results [20|, |2l| , previous works point to 
the presence of three regions: a region of extended eigen- 
states, a strictly localized region, with exponentially lo- 
calized eigenstates, and a mixed region, exhibiting both 
localized and extended features. The results of the litera- 
ture [201 suggest that the eigenstates decay algebraically 
in the mixed region. From the study of the IPR one 
can not distinguish between these two different types of 
localized states. 



IV. CONCLUSION 

In this paper we have studied the localization of eigen- 
vectors of random matrices through the calculation of the 
inverse participation ratio with the cavity method. We 
have found a self-consistent equation for the inverse par- 
ticipation ratio in the limit N -^ 00, which can be solved 
numerically through a population dynamics algorithm. 
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FIG. 6: Population dynamics results for the critical line sep- 
arating localized (L) from extended states (E) in the fully- 
connected Levy matrix. The results were obtained consider- 
ing e = 0.0001, 7 = 0.008 and a population of 10^ samples. 



Therefore, this approach contains no finite size effects in 
contrast with numerical diagonalization methods. The 
resultant equations for the inverse participation ratio are 
conjectured to be exact for Laplacian matrices on sparse 
random graphs and for fully connected Levy matrices in 
the limit iV — > oo. 

We have calculated the inverse participation ratio of 
Laplacian matrices on sparse random graphs. The spec- 
trum is characterized by a delocalized part centered 
around zero and a localized part in the edges of the spec- 
trum. The states corresponding to large eigenvalues are 
localized on single sites with large degrees in compari- 
son to the average connectivity [lal ■ Numerical diagonal- 
ization results for the inverse participation ratio appear 
to converge to the corresponding theoretical values for 
N -^ oo when there is no degeneracy in the eigenvalues. 

The matrix elements of random Levy matrices are 
drawn from a distribution with power-law tails character- 
ized by an exponent a. In previous works |20l.l2ll|. the lo- 
calization properties of Levy matrices have been studied 
using mainly diagonalization results. These results indi- 
cate that the eigenstates of Levy matrices undergo a tran- 
sition from a delocalized to a localized phase. However, 
large finite size effects are present in such calculations. 
Using the cavity approach we have determined a transi- 
tion line separating localized from delocalized states in 
the (a. A) plane by studying the behaviour of the inverse 
participation ratio in the limit A^ — > oo. Our results 
confirm the presence of a region with localized states for 
large eigenvalues, where the states are localized on pairs 
of very strongly interacting sites [20| ■ 

While we were writing this paper a prep rint appeared 
on the arxiv addressing similar issues |35| . The authors 
of this paper have determined the location of the An- 



derson transition in electronic systems on Bethe lattices 
using the cavity method, while we have focused on the 
localization properties of random matrix ensembles. 



Appendix A: The eigenvalue-dependent IPR and the 
Green function 

In this appendix we show how the eigenvalue- 
dependent IPR can be expressed in terms of Gf({X — ie) 
by means of eq. ^. By substituting eq. dl]), one can 
write down the following quantity 



gf{X) = \ime\Gf^{X-^e)\\ 



(Al) 



in the form 



N 



gfiX) = luney ^'^■'^^' 



-(A-AJ2+,2 



lim Df (e. A) 



(A2) 



where we have defined the non-diagonal contribution 

N 



-f(..^)-ET^^E '*''' 



M=l 



X— Xn + ie ^-^ X — Xu — ie 



(A3) 



Assuming there is a set of eigenvalues A = 
{Afcj, . . . , Afc^} such that A = A^. for any Afe. € A, we 
can rewrite Df [e, X) as follows 






I \2 



ii^l) 



X — Xu — ie 






+^E 



X — Xu + ie 



E 



,^A^-^^^+''.^.A^-^--'' 



(A4) 



In the absence of degenerate states in the spectrum, the 
set A is simply given by „4 = A^ , which reads 



u=jik 



J ^2 



I "' ^ X- Xu - ie 



+ ^E 



u^k 



X — Xu + ie 



E 



(V4 



J ^2 



A - A^ + ie ^ A - A^ - ie 



(A5) 



Thus we obtain Mm^^o Df {e, X) = and eq. (|A2p as- 
sumes the form 



N 



C(A)=^5](^;r<5(A-A, 
ti=i 



(A6) 



By summing the above equation over all the sites and 
dividing by Np{X) we obtain, in the limit N -^ oo, the 
identity that relates the eigenvalue-dependent IPR with 
G^(A — ie) (see eq. ©)■ The IPR associated to the state 



II is defined as Y^ = Eili(V'^)^- The identity §^ holds 
only in the absence of degenerate states. In the presence 
of degenerate states, the function &f (A) is given by eq. 
(jA6p plus a correction term that involves a sum over all 
the degenerate eigenvectors. 



[3l| . First, we assume that the functions Vj^\,{x) and 



,(i,fc). 



,W 



■ N 2 \^) factorize over the sites according to 'PVz(^) = 



tN 



,(0 



,(i,fe). 



tN 



n;;a.^]v:.(^.) and vz>{x) = n;;aAfe^)^;^ (^^)' ^^- 

spectively. Second, we assume that, in the limit N — >■ oo, 
the marginals on the cavity graphs fulfill Vj^^lxj) = 

When inserted in eqs. (|B2[) and (|B3|) . the above as- 
sumptions give rise to 



Appendix B: The cavity method 

Wc show in this appendix how to derive the cavity 
equations for the normalized complex function 'Pn,z{x) 
defined by eq. pT|) . We focus here on the ensemble 
of fully-connected Levy matrices in which Ja = for 
Vi. The values of the nondiagonal elements of J are 
i.i.d.r.v. drawn from the Levy distribution, defined in the 
subsection IIIIBI 



'PN,z{xk) ^ exp ( --^zxl 



y(l) 



Y[ / dxj V}^'., {xj ) exp {ixk Jkj Xj ) (B4) 



'PnA^i^) --exp {--zx,. 



X Yl dxj'P}!^^^{xj)exp{ixkJkjXj) . (B5) 



jedk\i 



1. Cavity equations 

The marginal at site k is defined as follows 



'PN,z{xk) = 



n '^•^j 

J&dk 



VnA^) : 



(Bl) 



in which dk denotes the set of indices in a row k for which 
Jij 7^ 0. Here dk is composed of a number of indices of 
0{N). 

Using eq. (|B1[) as a starting point, one can derive the 
following equations 



'PN,z{xk) 



n '^^^ 



K'li^) 



N 



X exp I zx/. + ixk 



/ , 'JkjXj 



(B2) 



oedk 



-P^i^k) 



n ^^j 



]edk\i 



-,{kd] 

N,z 



(X) 



N 



X exp --zxl + ixk ^ Jkj 



(B3) 



3edk\l 



where dk \ I denotes the set dk without site I, and the 



function V 



{ii,...,iM) 
N.z 



{x) is defined on the cavity graph 



g{ti.-,^M)_ The cavity graph t/(»i--*J'') is the subgraph 
of the original graph Q, in which the nodes ii, . . . , im and 
all their links with the other nodes have been removed. 
In order to close the system of eqs. (|B2p and (|B3p . 
we make two assumptions which have been used in the 
context of the cavity method for disordered systems 



The form of eqs. (|B4[) and (jBSp suggest that they can 
be solved in a self-consistent way through a Gaussian 
assumption for the functions T'j^{,{xi) and VN.,z{xi). The 
variances of the local marginals VN,z{xi) are the diagonal 
elements of the Green function, as one can note from eg. 
(|13p . Thus we make the following Gaussian ansatz [13[ 

for the cavity functions 7^Vz(^fc) 



T^N,zi^k) - ^1 ^ ^jv^(;) 



■ exp 



ixt 



2nG-,->{z) ' V 2G^*'^(-). 



(B6) 



where G/^i^ (z) are the diagonal elements of the Green 
function in which row / and column I have been removed. 
The substitution of the ansatz (jB6p in eqs. (jB4[) and (|B5|) 
leads to the following self-consistent system of equations 



G^^''\z) 



1 



^-9N,k(^) 



(B7) 
(B8) 

for z = 1, . . . , A^ and for all k <E di. The functions g%\{z) 
and hj^ (z) are defined as 



""■"^^^^^^^5^ 



mk 



W= E J^G^^'^i^) 



(B9) 

jed,\k 

h^[z)=Y^4G^^^{z). (BIO) 

The fixed-point solution of eqs. (|B7p and (jBSp allows one 
to determine the diagonal elements of the Green function 
for a single instance of J, which give access to the DOS 
and the IPR. 



2. The ensemble average 



where a'i is the variance of the weak matrix elements [32 



In this subsection we explain how one can perform the 
ensemble average and derive a self-consistent equation 
for W\^^{(jj) by employing a method introduced in [3j, 
|34| for a fully-connected Levy spin-glass. The method 
consists in the introduction of a small cutoff 7 that makes 
a distinction between small and large matrix elements. 
The global contribution of the small matrix elements is 
taken into account by means of the law of large numbers. 

We define the sets of indices that distinguish between 
weak and strong matrix elements in a certain row i ac- 
cording to 

0(7) = {] e Nn [i,iV]|(J,, < 7) A [j ^ z)} , 
Ui) = {j e Nn [i,iV]|(J,, > 7) A (j ^ z)} . 

These definitions allow us to rewrite g%\.{z) and hy/{z) 
as follows 



i6Ci(7) 



- 2^ JijGjj ' (z) 

jeCi(7) 






k e G(7) : 



ieCi(7)\fc 



7 C 



NX^) 



(^) 



E "^^j^j/ ' (^) 



E j^^G^r^ 



jeCi(7) 



iec,(7) 

E 

j6Ci(7)\'^ 



7-2 (^NXi) 



(z). 



In the limit iV — > 00, we can remove the k dependence 
from the sum over the weak matrix elements because it 
contains an infinite number of terms. Defining the joint 
distribution ^i {lo) of the real and imaginary parts of 
G^^\z) for a fixed j 



, N 
f|(^)(cj)= lim — V^[' 



-G 



N.U) 



(^) 



(Bll) 



one can apply the law of large numbers to the contri- 
bution coming from the weak matrix elements, leading 
to 



lim 



jeCiii) 



^ 4G;^-«(z) = aMdc.f^«H.., (B12) 



7 

N I dJPa{J)J^ 

-7 



2-a 



(B13) 



with the distribution Pa{J) defined by eqs. (|^ and ((^ . 

A self-consistent equation for Viz' (w) is derived by sub- 
stituting eq. (jB7P in (|BIip . One has to distinguish be- 
tween two cases. If j € 0(7): the function g}^' Az) is 
equal to h}^{z) for any index i, and as a consequence 

one obtains ^z {(^) = W\^Juf). If j € di.!)-, one 
can follow the discussion of [S^ in order to show that 
{li^\uj) = VFA,e(w). The distribution Wx,A^) fulfills the 
self-consistent equation 



Wx,.A^) - £ 



e '="' c;5; 



A;=0 



fc! 



]^dWiI^A,e,7(w;) 



1=1 



WdJipj^^{Ji) 



.1=1 



y. 5 \ to 



^7^ -J2l = lH{^l,Jl] 



(B14) 



in which (/(w)) = / (ia;iyA^e^-y(w)/(a;) and 

i/(a;, J) = J^w , 
2C„ 



0:7" 

_;,-^ |j|>7 

|J|<7 



PJni.J) 



The quantity Ca is defined by eq. ([28 



(B15) 
(B16) 

(B17) 
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